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Abstract 


We study a statistical ensemble of a single polymer with self gravitational interaction. 
This is a model of a gravitating string — the precursor of a black hole. We analyze 
averaged sizes by mean field approximations with an effective Hamiltonian a la Edwards 
with Newtonian potential as well as a contact repulsive interaction. We find that there 
exists a certain scaling region where the attractive and the repulsive forces balance out. 
The repulsive interaction pushes the critical gravitational coupling to a larger value, at 
which the size of a polymer becomes comparable to its Schwarzschild radius, and as a 
result the size of the corresponding black hole increases considerably. We show phase 
diagrams in various dimensions that clarify how the size changes as the strengths of 
repulsive and gravitational forces vary. 
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1 Introduction 


It has been long known that a typical configuration of a highly excited free fundamental string has 
a size R ~ where A^(S> 1) is the excited level and is is the string length scale. Since the mass 

of the level M states is proportional to y/M, which may be identified with a typical length of the 
string, the profile can be understood as that of a free random walk [T]0 The size of the configuration 
will decrease if the interaction between each part of a long string is turned on. At a sufficiently 
strong coupling, the size would be the string scale, and at that point Susskind conjectured that 
the collapsed configuration of a string may be identified with a small black hole of the same energy 
[3 ID- A natural question to be addressed is how the size of a long string changes as the string 
coupling constant varies. Horowitz and Polchinski [5] studied this problem by employing so-called 
thermal scalar theory, which provides the statistical nature of a long string near the Hagedorn 
temperature such as a spatial profile of the string^] Through a scaling argument, they estimated 
a typical size of a bound state wave function that describes a scalar field soared in a gravitational 
field of its own. 

We revisit this problem on gravitating strings with much more emphasis on a description based 
on ensembles of interacting random walks /polymersl§ Mathematical descriptions of random walks 
provide powerful tools to investigate statistical properties of polymers and so far lots of fruitful 
studies have been undertaken through methods such as the renormalization group [IIIIIIKIS]. The 
conformation of random walks is governed by two competing tendencies; the diffusion of monomers 
which is universal in systems at finite temperature and elasticity which is characteristic to long 
chains, both of which originate in the entropic nature of ensembles. The balance between these 
effects determines the characteristic features, such as the typical shape of polymers. Indeed, the 
balance between the diffusive and elastic forces results in the typical size Rq ~ A ^2 in an ideal 
random walk, where N is the number of monomers. A real polymer, unlike an ideal one, does not 
intersect with itself and is described by a self-avoiding walk (SAW). As is well known, the self- 
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avoiding effect makes the configuration puff up to a size of order N ^+2 in d spatial dimensions [IT] . 
This scaling is due to the entropic elasticity with the excluded volume effect taken into account. 

Though fundamental strings can be described by ideal walks, it has also been suspected that 
when it is compressed into a small region a self-avoiding property would emerge nonperturbatively 
m- Therefore it is intriguing to introduce a repulsive interaction in the description of strings as in 


^Manes has studied a Rutherford-type scattering for a free long string and provided a further evidence that the 
profile is indeed that of a free random walk [2]. 

^Properties of the thermal scalar are further investigated in [^, and recently more intensive studies have been 
carried out by Mertens et al. [3. There is another trial to evaluate it through the change of the density of states due 
to self interaction [8]. 

®We give a brief comment on introduction of a statistical ensemble. A string of a given mass y/M/is is specified 
by superposition of level M states. Since the number of states is enormous, a quantum average with respect to a 
randomly chosen state for macroscopic observables, such as a size, density distribution and so on, agrees with a 
statistical average over the level M states, which are dominated by configurations of random walks. Such a statistical 
average has also proven to be useful to derive a (pseudo-)thermal emission spectrum from long strings I10|. 
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real polymers and we will describe it as a type of excluded volume effect among monomersO The 
other types of self-interactions are also important. Such self-interactions include van der Waals 
interaction or electric repulsive interaction, for example, but our particular interest is gravity pro¬ 
duced by itself. Thus we shall investigate a long polymer with short-range repulsive and long-range 
attractive self-interactions as a model of self-gravitating stringsj§ The typical size of an interacting 
random walk is determined by interplay among the repulsive and attractive self-interactions in 
addition to the diffusive and elastic forces. Our main focus in the present paper is to see whether 
there exist any scaling regions in which gravity and the others are balanced, and furthermore to 
present an overview of the size behavior with respect to the change of the coupling constants of 
each self-interaction. 

This paper is organized as follows. In the following section, we shall argue the size estimation by 
using two methods, first by a variational method with harmonic potential, and second by a uniform 
expansion model. The results are organized into phase diagrams in various spatial dimensions. In 
section [3l we give a summary of the result. Details of the calculation in the uniform expansion 
model and a brief argument for van der Waals interaction are presented in the appendices. 

2 Self-avoiding random walk with long range attractive force 

In the random walk model, a long polymer is described by a chain of N numbers of monomers 
jointed freely with bonds. The mathematical description for the statistical property of a self¬ 
avoiding polymer with interaction is given by Edwards Hamiltonian [18] in d spatial dimensions, 

d / (9R \ ^ 

^ /o ( 2 . 1 ) 

Here the potential term consists of long-range Newton interaction in d (> 2) dimensions as well as 
a point-like repelling force expressed by the delta function, 

H(R(a),R(a')) = - , (2.2) 

where and u are dimensionless coupling constants. £ is the (Kuhn) length of the bond between 
monomers, which will be identified with the string scale ig- We are considering an ensemble of 
highly excited long strings of level Af, and the corresponding temperature /3~^ is on the order of 
the string scale 1/f's, but its explicit value is not relevant to our analysis. This sets the scale of 
the analysis, and the coupling constants are measured with respect to this scale. Since the length 
along a string is proportional to Igy/M, the excited level is related to the number of monomers as 
N oc VAf. The Hamiltonian (j2.Ih can be understood as a continuum version of a discrete model 

■^Real polymers may be under the circumstances with various external forces like an electric potential, or interac¬ 
tions with solvent. We shall not argue these external effects in the present paper. 

®The idea has already been mentioned in [5], and there have been a couple of works based on this picture [1511161117] . 
Our analysis partly overlaps with m, but we here pursue more variety of scaling behavior with repulsive interaction. 


3 





(a bead-spring model with interaction), 

d ^ 

E mn,R„'), (2.3) 

n=l n^n' 

where R„ is the position vector of n-th monomer. In the interaction term, n = n' terms are 
excluded since they correspond to the self-interaction of each monomeric In the continuum version 
(EU), this self-interaction of each single monomer should also be excluded and we can understand 
the coupling constants are suitably renormalized ones (with appropriate regularization such as 
point-splitting). 

What we want to understand is how the size varies as the couplings and u change. In 
particular, we shall explore if it exhibits any scaling behavior with respect to the number of the 
monomers N. Our primal goal is to evaluate the end-to-end radius squared, 

( r 2 ) =1 j VIl{N) (R(A^) - R(0))%-^-^ , Z = j , (2.4) 

where the end point R(iV) is integrated over while the other end is fixed due to the translational 
invariance. Hereafter, R(0) is fixed at the origin. 

Before starting the analysis, we roughly estimate the N dependence of the coupling constants 
with which the typical size of long polymers gets affected. As is well-known, with no interaction 
= n = 0, a typical size of free long polymers is given by Rq ~ (hereafter Rq always stands 
for this free walk size). When the coupling constants are turned on, the size will be altered but the 
N dependence remains the same until the couplings reach certain marginal values. Such marginal 
values, go and Uq, are given by the condition that the interaction terms are of 0(1). Since the 
distance of two different monomers scales as |R((t) — R(cj')| ~ the gravitational force becomes 
significant at go, where N‘^g‘^N ~ ~ 0(1) and the first comes from the number of pairs. This 

condition gives go ~ N~^ [5l[16]. On the other hand, the repulsive interaction is local. The volume 
occupied by the polymer is proportional to A" 2 , and then the chance of two different monomers 
coming to the same point is A“ 2 . Hence, the marginal coupling is determined by N‘^UoN~^ 0(1), 

d-4 

OT Uo ^ N 2 . This result suggests that for d > 4 the contact interactions become negligible, and 
this is consistent with the well-known fact that the probability of self-intersection of free random 
walks is negligible for d > 4. So far, we have investigated the marginal couplings at which each force 
takes effect against the entropic elasticity. It is possible that the gravitational and the repulsive 
forces balance out to sustain a configuration. This is realized when g'oN ~ ~ ~ (the 

®Precisely speaking, a free part should have been written as ^ IH-n — Rn-i|^ with a fundamental bond 

length b and a total number of monomers N. The continuum limit is taken by 6 0 and iV —>■ 00 with the total 
“diffusion time” T = Nb^ fixed, and the parametrization on the polymer is t = nb^ (0 < t < T). (This is the 
continuum limit for the diffusion equation that a random walk probability distribution satisfies.) Kuhn length £ and 
an effective number of monomers N are introduced so that the size, the total mass and length be £N'^, N/£ and N£ 
respectively; namely a polymer is effectively regarded as N number of monomers (of unit mass) joined by bonds of 
length £. T is identified with Nl'^, and a continuum dimensionless variable a = tjR (0<(j< N') is introduced to 
represent a position on the polymer. 
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prime is put to indicate that a different type of equilibrium is realized), or Ng'^/u'^ ~ 0{1). We 
will see that the following calculation reproduces this condition. 


2.1 Evaluation of size by variational method with harmonic potential 


We shall evaluate the averaged size-squared (R-^) by the variational principle. The trial Hamiltonian 
is chosen to be a harmonic action, 

2 


‘'"(f) 

where g is a dimensionless variation parameter^ The free energy, 

pN pN j 2 pN 

f3{H-Ho)=J^ da da'V{B.{a)M^'))-i^ 
satisfies the following inequality, 

/?F</ 3Fo(g) + /?(R-Ro)o , 


(2.5) 


( 2 . 6 ) 


(2.7) 


where (• • •)g is the expectation value with respect to the trial Hamiltonian (5Hq and /3To is the 
corresponding free energy. The variation parameter q will be chosen so that it minimizes the right 
hand side of the inequality, and we shall estimate averaged sizes by using the optimized parameter. 
The trial Hamiltonian provides the propagator 


G{a,a') = 


qd 


27rt'2 sinhg|cr — a'\ 


exp 


qd [(R((t)^ -|- R(it')^) cosh( 7 |cT — a'\ — 2R((t) • R((t') 


2£2 sinhg'Icr — a'\ 


( 2 . 8 ) 


The expectation value of a function of k different points R(cri) with respect to PHq is calculated 
by use of G{a, a') as 


(0(R(a),--- ,R(cTfc)))g = ^ 


e-dHo 
k+l 


Y ffl ['^R(^i)G(ai.i,a,}]0(R(a),--- ,R(<Tk)) , (2.9) 

° i=i 


with (To = 0, (Tfc+i = N, and the partition function is 

Zo = y PR(iV) G(0, N) = cosh"i qN. 

The expectation value of the size-squared with respect to /3Ido is 

If ^ 

(r2)^ =_ / VR{N)R{N)^G{0,N) = — tanhgiV. 


( 2 . 10 ) 


( 2 . 11 ) 


^This corresponds to the one used in m cLS ^ours — 
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In the following, we shall find optimized parameters, which we call qo, that minimize the free 
energy bound (I2.7jl . They will be given in terms of u and g, and determine the size behavior in the 
space of the coupling constants. We claim that the mean radius of configurations which we 

shall denote R^, is approximated by using the optimized parameter as 

~ (R^) = — tanhgo^^- (2.12) 

^ qo 

From this expression, one can see that the size stays to be the free walk one Rq for qqN <C 1 and 
it changes from qqN ^ 0{1) as 


R ~ < 


iVN qoN^l 
£ 


(2.13) 


- y/% 


qoN > 0(1) 


Now we evaluate the term {f3{H — Hq))q in the right hand side of (|2.7p . For the quadratic 
potential part, we find 

j2 rN 

^ . II III iti'n I 1111 I \j 

(2.14) 


2P 


j /-or ^ 2 \ qdiPPo) qdN 
da (R(ct) = 2 ^^ = ^tanhgiV. 


As for the original potential term, we first note that the part including the Newton potential can 
be rewritten in terms of an exponential, 

d-2 

r{^)Jo \ /c’ ^ ^ 


1 


|R((t) — R((t' 


/\\d-2 


0 


0 


where e > 0 is a regulator. The integrals with respect to R(cr) are just Gaussian. Therefore a 
straightforward calculation gives the expectation valu^ 

1 \ 1 


R((T)-R(a')|‘'-Vo r(|) \2£^F,{a,a'-,q) 


qd 


d-2 

2 


(2.16) 


where 

f sinhg'iTCOsh(7(A'— (t) + sinh(7(T'coshg(A^ — cr') — 2sinhgcJCOshg(A'— a') 

For the repulsive interaction part, we rewrite the expectation value of the delta function in the 
integral form. The integrals are again Gaussian, and we easily find 

d‘^k 


(5('')(R(u)-R(u')))^ = I 


(2vr) 


,*k-[R(<T)-R(cr')]\ _ 


qd 


0 \2TTi‘^F2{a,a']q) 


where 


F 2 {cr, a';q) 


sinh qa sinh q{a' — a) cosh q{N — a') sinh qa' / sinh qa 


sinh qa' 


+ 


cosh qN 


1 - 


sinh qa' 


, (2.18) 


(2.19) 


® In what follows, we take cr < ct' by writing the integral as da da' = 2 da' da, since the integrand 
is symmetric under the interchange of a and a'. Then the results do not appear in a symmetric way for exchanging 
a and a'. 
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Combining these into (j2.7p . we find the inequality 


/3F < ^ In (coshgA^) — tanhgA^ 


-2 


rN nc 

/ ^^7 

JO JO 


da 


9 


ni) 


qd 


2Fi(cr, cr';g) 


d-2 

2 


— W 


qd 


27rF2{a, a'-,q) 


( 2 . 20 ) 


We need to tune q to find the most strict bound for the free energy. It is, however, difficult to 
analytically evaluate a and a' integrals, since Fi{a,a';q) and F 2 {a,a']q) are complicated functions 
of a and a'. We instead use an approximation to proceed. We take all the exponential quantities, 
namely and to be small and negligiblel^l This will 

not hold near the ends and the point at which cj = cr' in the integration region. These points are in 
reality separated by a fundamental bond length, and the approximation will be valid elsewhere as 
long as q is not small enough. With this approximation, the free energy bound is simplified to be 


/3F < 





( 2 . 21 ) 


With unimportant positive numerical factors omitted, the extremal condition becomes 


l-Ng^q'/ +Nuq"/ =0. 


( 2 . 22 ) 


In the following, we examine (j2.21j) and (j2.22p to find an optimal value go- We start with a pure 
gravitational case {g > 0,u = 0) to see the method rederives the already known behavior in [5l 116]. 
Next, a pure repulsive case (u > 0,g = 0) is argued and a limitation of the current variational 
method is presented. Finally, we consider a generic case {g, u > 0) in various dimensions and 
discuss that there appear two different size scalings with respect to the magnitude of u. 

We first consider the case with no repulsive force, namely u = 0. This corresponds to the 
original self-gravitating fundamental string studied in um- The results depend on the spatial 
dimension d. For d other than 4 (for which a separate argument will be given shortly), the extremal 
value is go — (a numerical coefficient is of no importance). In d < 4, the right hand 

side of ()2.2ip is a convex function of g when g^ > 0, and then go indeed provides the optimal value 
(an absolute minimum), and the marginal coupling go is obtained form the relation q^N ~ 0(1) as 

d—6 

go — N^~ which is consistent with the previous estimation at the beginning of this section. The 
size behavior is easily read off from 


R ~ 


®With this approximation the propagator 


Gofer, cr') ~ 


dq \ ' 
ttF J 


exp 


iy/N 9 <90 

i[g‘^N)'^ g>go 

becomes that of the ground state approximation, 


^ (R(cr)^ -h R(cr')^) - y |cr - cr' 


(2.23) 


Thus, this can be viewed as the application of the ground state approximation to the evaluation of the right hand 
side of (|2.20l) . which is valid if the dimensionless level separation q is not much small. 
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The size would change as the coupling grows and becomes comparable with the Schwarzschild 
radius Rg — at the critical coupling gc ^ N 2 ^ where we have identified ^ with the 

string fundamental length 

In d = 4, the free energy satisfies the inequality fiF < qN{l — g‘^N'j. Note that the marginal 
and the critical couplings are the same order, namely go ^ Oc ^ N~ 2 . When g ^ go, the optimal 
value is qq = 0, which gives the free random walk behavior. On the other hand, when g > go, we 
find go ^ oo and i? —)• 0 as is obvious from (12.131) . This means that the polymer suddenly collapses 
at go- In the case d > 4, the marginal coupling go is larger than gc, and the system is suspected 
to exhibit a hysteresis [S]. Since the right hand side of (I2.2ip is a concave function, the extremal 
value go — leads to a local maximum, and the correct optimal value of g is go = 0 for 

= 0 and go = oo for > 0. This implies that the configuration completely collapses from a 
free walk size Rq with an arbitrary small value of g. The behavior is more exotic than discussed 
in [S] (this is also mentioned in [IB]). In this pure attractive case (especially 2 < d < 4), we have 
rederived the results obtained by Horowitz-Polchinski [5] and Khuri [El- 

Next, we examine the case with a repulsive interaction and take u ^ 0. To begin with, we 
give a comment on the case with a pure repulsive force m > 0, g^ = 0. This should correspond 
to a SAW without the attractive interaction, and one may expect to obtain the famous result by 

_ 3 _ 

Flory [m, no Nd +2 (1 < d < 4). From (I2.2ip one can immediately see that the only solution is 
go = 0 when g^ = 0 since all the coefficients are positive 0 This provides the free walk behavior 
of the size Rq and fails to reproduce Flory’s result. Indeed, as long as we consider a confining 
harmonic potential, the size of the walk should decrease compared to the free walk result for any 
choice of g. One sees go = 0 stands for no harmonic potential, so the configuration expands as a 
free walk, but it cannot spread larger than that sizelHl Based on this observation, we conclude that 
the variational calculation with a harmonic potential is not suitable for describing random walks 
that might expand larger than the free walk size Rq. From (|2.13p . one sees that the size reaches to 
the fundamental length scale at go ~ 0 ( 1)5 and the description based on this effective Hamiltonian 
may not be valid after that point, providing the upper bound on go. Thus, the validity of this 
variational method will be as follows: if the maximal size of a configuration is known, a priori, to 
be Rq, we may trust the method for go < 0(1), but if it is possible that a configuration can expand 
larger than Rq, it is safe to take go to be N~^ ^ Qo ^ 0(1). At any rate, this variational method 
based on the harmonic potential will serve a reasonable size evaluation if the size is smaller than 
or equal to Rq. 

_ d-4 

Now we turn on the gravitational force, g^ > 0. The only negative term in (12.220 . —Ng'^q^^ , 

d-2 

will balance out with the dominant one between the positive terms, 1 and NuqQ ^ (or both if they 
are comparable). If the gravitational force term balances with the first term, 1, the solution is 

^'^Note that go has to be non-negative otherwise the expectation values are not well-defined. 

^^This argument can also be applied to the situation with sign flipped, —>■ —g^- Together with a repulsive 
force u > 0, this describes a single polyelectrolytes chain in an ideal situation where no screening effect from solvent 
takes place. By the Flory’s scaling argument the size behavior is evaluated to be i? ~ N ^ [19], and furthermore by 
employing a renormalization group argument R ~ N 4-2 [lailSj; they are much larger than the free walks, but our 
calculation merely gives a free walk behavior. 













go ^ as in the previous pure gravity case (we give a separate consideration for the 

d-2 

d = 4 case later). This solution is valid as long as NuQq^ <C 1 for a given value of go- This 
condition is rephrased in terms of a condition for g, 

g<go (2 < d < 4), (2.24) 

where 

go ~ tt2(d-2) pj~~d=2 . (2.25) 


In d > 4, the condition for this go to be valid is g > go, but the configuration with pure gravity 
(namely, without the repulsive force) exhibits a strange behavior as we have seen, and the variational 
method itself does not seem valid. On the other hand, if the third term in (|2.22l) is dominant over 
the first term, the gravitational and repulsive forces balance, and the solution is go — ^ ju. The 

d-2 

consistency condition, Nuq^^ ^ 1, implies that g < go for d > 2. Thus, go is a marginal coupling 
at which the solution of the stationary condition switches from go — (g'^A^^) to go — g'^ju. At 
g = go, all the terms in (I2.22p are of the same order, and these two expressions of go agree. The 
observation so far can be demonstrated by use of an explicit solution of go in d = 3, 


do = 


2Nu 


( - 1 + ^yl + 4N^g^u) 


(2.26) 


In d = 3, go c; Thus, g < go (in terms of N dependence) implies that N'^g^u <C I, and 

in this region go is reduced to go — g^N"^- On the other hand, when g > go, we have go — g'^ju. 
Around g go {N'^g'^u ~ 0(1)), the entropic elastic, the gravitational, and the repulsive forces 
balance out, and go ^ which can be expressed as g^A"^ or g'^/u by use of u ~ {gN)~‘^. 

We move on to the separate case in d = 4, where go ~ N~ 2 . We find 


do 


g^N-i 

uN 


(2.27) 


Note that g‘^N has to be greater than or equal to 1 since go must be non-negative. For g < go, namely 
g‘^N < 1, g^ term in (12.221) is g independent and subleading. We then come back to (I2.2ip and find 
that go = 0. For g > go, go becomes g^/w, and for g ~ go (g^A ~ 0(1)), do — {uN)~^ ~ g‘^/u. 
Thus, in d = 4, the optimal value of go changes from go = 0 to go = g^/n at g = go- 

Let us estimate the size scaling in the case of go = g^/tt in arbitrary d > 2. The mean square 
size is given by (I2.12p . 




i'^u 

~ —— tanh 
9^ 





fN 


g^Nju ^ 1 
g^N/u ^ 1 


(2.28) 


Thus, when the repulsive force is much stronger than the attractive force, u g^A, the size is 
given by just the free random walk one, but does not expand larger than that. The feature that 
the size does not exceed the free walk size would be a limitation of this approximation as we have 
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discussed in the pure repulsive case. In the opposite limit, the behavior is affected by N dependence 
of the interactions. Here we consider the N dependence of u to be fixed, and observe the scaling of 
the size by changing g^. When qqN ~ 0(1), the size starts to change and this defines the marginal 
value 


9 


/ 

O 





(2.29) 


which agrees with the naive scaling argument at the beginning of this section. As the coupling 
grows, the size decreases, and it eventually coincides with the Schwarzschild radius of a black hole 
of the same total mass; Rs — in d dimensions. The value of the critical coupling at 

which this transition takes place is 


9 


/ 

C 


d-2 1 


CliU N d , 


(2.30) 


and the size of the corresponding black hole is 

RcC=t£{uN)^ . 


(2.31) 


Let us summarize the size scaling in this generic case. The size depends on the A^-dependence of 
the coupling constants g and u. Here, we observe the change of the size by tuning the gravitational 
coupling g from o' = 0 to larger values for a given value of u, until the size of the configuration 
becomes the Schwarzschild radius of a black hole of the same mass. For a small value of g' < ^o) 
the first two terms in (j2.22p balance, and the size scaling is the one without a repulsive force given 
in (j2.23p . Thus, for sufficiently small couplings g < go,9o, the size is the free walk one Rq. If the 

d-4 

repulsive force coupling is larger than the marginal value Uq ^ N 2 ^ the configuration is expected 
to be expanded larger than Rq, but the variational calculation is not capable of realizing such an 
expanded conhguration as argued. In order to obtain a consistent picture, we take u < Uq with 
which the size is to be Rq for a very small g < go, 9o, and leave the analysis for u > Uq with small g 
to the next subsection where another approximation method is introduced. For larger values of g, 
the behavior depends on the magnitude relation among go, go, and g'^. It is not difficult to check 
that g'^ < (jo u < Uo- Thus, in this region of u, the configuration immediately starts to shrink 
as i? ~ l^fujg once the repulsive force participates in the balance at g = go- go < (jo, the size 
shrinks as i? ~ i{g‘^N)'^ from g = go and changes its behavior to i? ~ ly/ujg at g = go, and the 
conhguration will eventually be covered by the Schwarzschild radius Rc at g = g^- However, if u 
is too small, the conhguration becomes a black hole before the repulsive force becomes important. 
It happens when gc < Qp , which leads to tt < N~^, and the conhguration collapses to a black hole 
of radius i at g = 5cO If (jo < 9o, the size starts to change as i? ~ i^fujg at g = go from the 
beginning, and it becomes a black hole at g = g'^- 

To illustrate a typical size behavior, we pick up a couple of examples in d = 3. It u takes 
the marginal value Uo ~ N~^, the repulsive force pushes the critical coupling to a larger value 

^^From (12.131) . one finds that the size becomes comparable to a fundamental scale £ at go ~ 0(1). If u < N~^, the 
u term in (12.2211 remains subleading for go < 0(1), and then negligible in the whole process. 
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g'^ ^ N 12 than that of non repulsive case, ^ N 2 ^ and the black hole swells {Rc — iNe). 
As we will see in Sec. mi real polymers have u ~ 0(1); which corresponds to a well-known 
phenomenological size scaling by Flory. If u ~ 0(1) in d = 3, the critical coupling and the size 
at the corresponding point are even larger, g'^ ~ N~3 and Rc ~ £N3 respectively^^ Note that 
from (j2.3ip one can see that Rc corresponds to a close-packed configuration of N monomers in d 
dimensions for rt ~ 0(1). In d = 4, ~ ^o, and one can see that the size suddenly changes from Rq 

to d? ~ i^fujg at g = (7o ~ N~^. Especially, if u < N~^, the configuration collapses in to a small 
black hole of size 1. More detailed observation on the size and the values of the critical couplings 
are presented in Sec. 12.31 where phase diagrams are drawn. 

So far we have obtained the size scalings which are reliable for the configuration which shrinks 
due to the attractive force. In order to investigate configuration larger than Rq, we need to invoke 
some other methods. In the following subsection, we shall consider a different approximation scheme 
that turns out to work for expanded configurations. 


2.2 Evaluation of size by uniform expansion model 


In this subsection, we try another approximation to evaluate the change of the size of self-gravitating 
polymers. This method is called the uniform expansion model (UEM) [13], in which the fundamen¬ 
tal length of the bond i is renormalized as i' = ai {a > 0) but the configuration itself is assumed 
to remain that of the free random walk. 

We consider an effective free Edwards Hamiltonian with the bond length T = ai, 


pH' =- 


d 


rN 


dal 


(S)' 


2a^i^ 

The Green function corresponding to this Gaussian action is 


G\a,a') = 


d 


^TxiPaPla — a'\ 


exp 


d 


2i‘^a‘^\a — a 


-(R(u)-R(uO)' 


(2.32) 


(2.33) 


and the expectation value with respect to PH', {■ ■ ■)', is calculated with this Gaussian propagator 
in the same way as (12.9p . where the partition function is unity. 

The end-to-end distance squared, evaluated with respect to the original action including the 
interactions m, is rewritten as an expectation value with respect to PH' as 


/(R(iV)-R(0))^e-^'^ _ ((I^(^)-R-(0))^e 

In UEM, we assume that P(H — H') can be treated as a perturbation, and then 


(r2> = 


(2.34) 


^ ' {l-P{H-H'))' 

noted a couple of times so far, when u > Uo, the variational method does not give a reasonable size scaling 
as long as g is small, but the critical coupling and size here are indeed valid as they agree with the result by another 
method we are about to introduce. 
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~ ((R(iV) - R(0))2)' (1 + {fi{H - H'))') - ((R(iV) - R(0))2/3(R - H'))' , (2.35) 

to the first order. One can easily see ((R(A^) — R(0))^)^ = , since it is the radius-squared of 

a free random walk with . The other terms are spelled out as 




36) 


and 


((R(1V)-R(0))2/3(R-R'))' = (R-(^)-R(0))'y 

rN pN 

+ / da da'{V{a,a'){Ii{N)-KiO))^)' . (2.37) 

Jo Jo 

We evaluate these expectation values in Appendix [X] in which the details are presented. The result 
of the expectation value of the size squared is summarized as 


(r2) =Ne'^ + Na^f{l - a^) + 2ufa 


2g^f‘aJ 

ni) 


/ 1 \ po 

Jo 


/ / \ 

da (a — a) 2 


(2.38) 


The integrals should be understood as the double sum, Yln'= 2 Yln=i^ much easier to 

evaluate them by using continuum variables a and a'. The result becomes 


(r2> = \Na^ + Aro2(l - a^) + CiuN^a^-^ - C2g^N^a 


4-d 


where 


Cl = 2 


d\2T{^) 


27rJ r(^)’ 


C2 = 


d — 2fd 


d 


■ 2r(^) 


(2.39) 


(2.40) 


are positive constants. They are divergent at some even d, but this comes from the point a = a' in 
the integral and the divergence is an artifact due to using a continuum variable. In evaluating the 
size, what we need is not the precise values of Ci and C 2 , but the fact that they are ^independent 
positive constants. We shall regard them simply as positive constants and omit them from the 
expression in any case. 

The consistency condition for UEM is that the size should be given by the free walk of the 
fundamental bond size i', (R^) = namely the parameter a is chosen so that the 

second, third and forth terms in ()2.39p cancel out. The consistency condition becomes 


a " - a '^+2 + = 0 , 


2„2 i 


6—«i 


(2.41) 


where we have dropped unimportant numerical coefficients. Note that if = u = 0, we have 
essentially a unique solution a = 1 and we recover the free walk result. In the following, we shall 
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find the solution of a in various situations. First, a pure repulsive case (u > 0, (7 = 0) is discussed 
to be capable of realizing expanded configurations larger than Rq] especially Flory’s exponent for 
real polymers is successfully reproduced. We briefly mention that UEM is not suitable for a pure 
gravitational case (5 > 0 , tt = 0 ), and then move on to the discussion of the generic case ( 5 , u > 0) 
where a size scaling similar to the one in the variational method is found. 

We start with the case with no gravitational interaction, namely a real polymer ( 5 ^ = 0, « > 0). 
The consistency condition is = 1 + ua~^N^~. If the repulsive force is very weak {uN^~ <C 1), 
we may take the leading order in u, and find ~ + uN ^~). Since S> 1, as we have 

discussed at the beginning of this section, the size is sensitive to the value of u only for d < 4. On 

4-d 

the other hand, for uN 2 3> 1, the condition leads to 

(2.42) 

and thus we find 

(2.43) 

d—A 

At the marginal coupling Uq ~ N^~, the size becomes Rq. When u ~ 0(1), we have the relation 

3 _ 

R ~ iN which is known as Flory’s exponent for real polymers [TI]. As just mentioned, this 

result is valid for d < 4, and we have witnessed that UEM is capable of evaluating a scaling size 

which is larger than that of free walks. 

Next, we consider the case with pure attractive force {u = 0, > 0). The self-consistency 

relation is a = 1 — g a N^~, For g N^~ <C 1, we observe a perturbative correction to the 
size, and find R^ ~ N£‘^(^l — g‘^N^~ ). Erom this relation, the marginal coupling go again appears 
to be go ~ N^~ [5]. Since a has to be positive, there does not exist a solution at strong coupling 
g‘^N^~ ^ 1 (namely g > go)- Thus, UEM fails to reproduce the scaling by Horowitz and Polchinski 
[5]. If we switch to a repulsive long range force, g'^ < 0, we find 

~ {-g‘^)-dN^ , ^ . (2.44) 

For \g‘^\ ~ 0(1), this result agrees with Flory type of scaling argument for a polyelectrolyte (a 
single charged polymer) [19], though it is not exactly the same as that of renormalization group 
analysis R ~ iNd^ [laiii]. 

Finally, we consider a general case {g'^ > 0,u > 0). The number of terms in the consistency 
condition (j2.4ip that come into balance at the leading order in N would vary. As in the analysis 
of the variational method, we take u as a given value and observe the size change as a function 
of g. Let us first consider the situation in which the repulsive force is effective, u > Uo- When g 
is small enough, a and the size R are given by (12.421) and (I2.43p respectively. The size starts to 
change when g term in (I2.4ip is comparable to the second and the third term as 

a^+‘^ ~ uN^ ~ g‘^a^N^. (2.45) 

d 3 

This leads to a marginal coupling g" — u'^id+ 2 ) jy d +2 ^t which a starts to decrease; for larger values 
of 5 , it is easy to see that the last two terms in (I2.4ip (or (I2.45h ) are dominant. This determines a 


13 















as 



which leads to the mean size squared 


(2.46) 


^ £ 2 ^ 2 ^ _ _ ^ 2 . 47 ) 
Note that this takes the same form as in the harmonic potential analysis (I2.28p . Thus, when 

1 3 

u > Uo, the size begins with an expanded size R ~ 1^^+^ N for small g, and as increasing g the 
configuration starts to shrink as R i^pujg at g = g", and eventually collapses to a black hole 
of size Rc at g = g'^ where the critical size Rc and the coupling g'^ are defined in (I2.3ip and (|2.30l) 
respectively. (This analysis provides u > Uq part of Fig. [1] and Fig. [3] given in the next section). 
Note that if u > the critical coupling g'^ is larger than the marginal coupling g" which implies 

that the configuration is covered by the horizon even before its size starts to shrink. 

Next, we consider the opposite situation, u < Uq- For a small value of g, a'^ and terms 
in (I2.4ip give a dominant balance solution a = 1. Thus the size starts with Rq. As g increases, 
the term g'^a?N~^ eventually becomes 0 ( 1 ) at g = go, and the solution of a starts to change. 

4 — d 

Since uN^~ <C 0(1), the repulsive force does not participate in the balance condition unless N 
dependence of a changes as a <C 1. However, before the repulsive force becomes effective, the size 
behavior is the same as the pure attractive force case ( 5 ^ > 0 , u = 0 ) ~ — g‘^N ^~) 

which has already been analyzed, and UEM fails to provide a consistent size decreasing behavior 
as discussed there. Therefore, we conclude that the UEM analysis is not reliable when u < Uq, and 
we use the result of the variational method in this region instead. 


2.3 Phase diagrams 

In this subsection, we organize the size behavior analyzed so far into “phase diagrams” in various 
dimensions. As discussed in the introduction, the shape of polymers is determined by the confronta¬ 
tion among four types of self-interactions; the entropic diffusive and elastic forces, the repulsive 
interaction, and Newton gravity. Two of them, the diffusive and the repulsive forces, make the 
size of configuration larger, while the other two, the elastic and Newton forces, are to compress 
the configuration. Thus, typically, there are four types of configuration that are characterized by 
the balance between two out of these four interactions; one from the former expanding forces and 
another from the latter attracting forces. These regions have different size scalings that depend on 
g, u, and N, and are separated by boundary lines that are determined by the conditions which two 
forces coming to balance. We call this diagram simply a “phase diagram.” The boundary lines are 
parametrized by various marginal couplings^ Uq ~ N^~ (the repulsive force becomes effective 
against entropic elasticity), go ~ N~^ (gravity starts to balance with entropic diffusive force), 

d—4 1 

§0 — u‘^(d-2) N~~^ (gravity, the repulsive, and entropic forces are in equilibrium), and g” (gravity, 

Another marginal coupling g'a — \/^, at which gravity starts to change the size against the repulsive force, is 
also important but does not appear as a boundary line. 
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the repulsive, and entropic forces come to balance in an expanded configuration). There also exists 
a special region in which the whole conhguration goes inside the Schwarzschild radius of a black 
hole of the same mass; we call it a “black hole” phase. The boundary of this domain is given by 
either of the critical couplings, Qc ~ N~^ (collapse against the entropic force) or ~ N~d 
(collapse against the repulsive force). 

Fig. □ shows the phase diagram for 2 < d < 4. The horizontal and the vertical axes are log^y g 
and log^ u respectively, and N independent factors are neglected. There are five domains that 


log AT W log AT i?c 



Figure 1: The phase diagram in 2 < d < 4. We set (. = 1. The right vertical axis shows the size of 
a corresponding black hole for a given value of u. 

we have just discussed, and we look at each of them in more detail in order. On the left-bottom 
corner, for small u and g, there is a “free polymer” phase whose size is given by Rq ~ £y/N in 
any dimension. In this domain, the diffusive and the elastic entropic forces balance out, and the 
size is stable against the change of u and g until they reach certain marginal values. Going up 
vertically, the configuration starts to expand at Uq, and we come into the region that we call “puff- 
up configurations,” where the repulsive force and the entropic elasticity balance out. The boundary 
is given by Uq- The size is given by i? ~ which is stable against the change of g but 

changes as u varies. We return to the free polymer domain and go to right as g gets larger (in the 
u < Uo region). It is not difficult to check that < do for ti < tio in 2 < d < 4, and Newton force 
comes to balance with the entropic diffusive force at go- On the right of this boundary, we are in 
a phase called “change I” in which the size decreases as g increases, as d? ~ £{g"‘N)d-'t ^ but is not 
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sensitive to the change of u. This is equivalent to the behavior of free self-gravitating polymers (no 
repulsive force), and it may collapse to a black hole at the critical coupling Qc- It indeed happens 
if u < N~^. When u > N~^, Qq is smaller than gc, and Newton force and the repulsive force will 
balance out before it becomes a black hole. Thus, after crossing a border line parametrized by a 
marginal coupling 5 ^,, the different size scaling R ~ i^fujg applies, and we call this region “change 
II.” Note that g'^ < go and the size immediately starts to change after crossing the boundary line. 
One can also come into change II from a puff-up configuration region, for a given u > Uq, by 
increasing g larger than g". Thus, the boundary between change I and change II is given by go, 
while the one between Puff-up and change II is by g”. If we further increase g in change II region, 
the configuration will be smaller than the size of the horizon at the certain critical coupling g'^. 

These phase boundaries are shown in Fig. [1] The green horizontal segment is given by Uo which 
separates a free polymer and a puff-up regions. Three blue lines (blue, light blue and cyan) denotes 
the point at which Newton gravity participates in force balance. They are parametrized by go, go, 
and g'o- The red lines consist of gc and g'^ at which the size of the horizon catches up the size of 
a polymer, and the configuration may collapse into a black hole. It is easy to see that g'^ is more 
steep than g” in general, and at log^v <7 = and log^yu = d — 1, 5 ' becomes smaller than g”. 
After this point, the configuration may become a black hole even no interaction work to shrink 
its size. On the right of the figure, the size of a corresponding black hole Rc for a given u is also 
shown. For u < N~^, the corresponding point is given by gc ~ N ~2 and Rc — i. If the repulsive 
force is significant in balance conditions, the size of corresponding black holes gets enhanced, as 
logj^{Rc/i) ~ ^ (-ggg (|‘2.3ip h and black holes swell in general. 

Fig. [2] shows the size change with respect to logjy g for various values of u in d = 3. This figure 
is schematic and the scalings are adjusted to draw the diagram. We pick up four typical values of u; 

3 3 

u ~ N ~2 (the repulsive force does not work throughout), u ~ N~i (the configuration experiences 
both change I and II), u AtO 

(a real polymer whose largest size is given by Flory’s scaling), 
and u ~ N'^ (the horizon size catches up before the interaction starts to work). The red line here 
shows the size of a horizon for the same mass. The crossing points with each line determine the 
corresponding points. 

Next, we quickly examine the cases in d = 4 and d > 4. From Fig. (TJ one can see the point at 
which four different regions meet (at (log^y g, logjy u) = (^^, ^^)) moves to a right-up direction as 
d comes closer to 4. At d = 4, some boundary lines merge and we obtain the left diagram in Fig. [3l 
The change I region disappears, and there are four regions left. For n < Uq ~ N^, the size is given 
by Rq in a weak coupling region. As g gets larger, it may collapse into a black hole if u < N~^ or 
may start to change its size as i? ~ iy/ujg. A crucial difference from the situation in 2 < d < 4 is 
that for u < Uo, hy crossing a vertical boundary line (denoted by thick light blue and red lines), 
the size of a configuration jumps; it suddenly collapse to a black hole of size i u < N~^, or it 
quickly shrinks down to the size R determined by the value of u for N~^ < u < N^. For u > Uo, 
the repulsive force is significant, and the behavior is analogous to the case in 2 < d < 4, and the 
size smoothly starts to change by crossing the boundary line g”. In d > 4, the diagram is given 
by the right in Fig. [3l It is similar to the situation in d = 4, but in this case the size scaling for 
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logiv^ 


Figure 2: The change of the size for various values of u in h = 3. £ is set to be 1. The red 
line represents the horizon size at g. On the horizontal axis, we show which marginal or critical 
couplings are relevant for each numbers. Note that go, g", and g'^ are u dependent. 

u < N~^ is ambiguous. As discussed in Sec. 12.11 the size scaling shows an anomalous behavior if 
the repulsive force is turned off. In u < N~^, the repulsive force is not effective and the situation 
is similar to that pure gravitating case. The other part is close to the d = 4 diagram. 

3 Conclusion 


In this paper, we have evaluated the averaged size of a long polymer which interacts with itself 
through Newton gravity and also has self-avoiding property that is implemented by a contact 
repulsive interaction, in spatial d dimensions. The mathematical description for the statistical 
property of polymers is given by a self-avoiding random walk with self-interactions, and is analyzed 
by an effective Hamiltonian a la Edwards. We have evaluated the expectation value of the end- 
to-end distance squared with respect to this Hamiltonian by employing two approximations; a 
variational analysis by a trial Hamiltonian with a harmonic potential and a uniform expansion 
model. The size of the polymer changes as the strength of the interactions varies. It has been 
found that the variational calculation gives a reasonable result when the size gets smaller than that 
of the free random walk Rq ~ IN 2 , in which the attractive force overcomes the repulsive or entropic 
diffusive force. This method, however, fails to explain the case where the configuration expands; 
even if the repulsive force becomes dominant, the size remains as Rq. This may be a limitation 
of the approximation due to using a harmonic potential to describe a dominant configuration. On 
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Figure 3: (Left) The phase diagram in d = 4. The vertical thick lines (with two colors) denotes the 
point where the size of configuration jumps. (Right) The phase diagram in d > 4. Below u ~ 
the behavior of the size is rather obscure. In both diagrams, I is set to be 1. 


the other hand, the uniform expansion model, where the configuration is fixed to remain the free 
walk and the change of the size is encoded in a renormalization of the length of the bonds between 
monomers, turns out to be capable of explaining reasonable changes of the size in the expansion 
case, where the repulsive force is sufficiently strong. As the size decreases, this method does not 
work at strong coupling for the pure attractive case, but it provides a consistent behavior when 
Newton and the repulsive forces balance out. In the most compelling case where the two forces are 
in balance, the size is found to be 

= (3. 

where g and u stand for the coupling constants of Newton and the repulsive forces respectively, and 
this scaling is valid for both contraction and expansion cases. The size scaling varies as which two 
of four competing self-interactions are in balance, and there in general exist four scaling regions 
depending on the values of g and u. Together with a situation in which the whole configuration 
is covered by the Schwarzschild radius of a corresponding black hole, unified picture of the size 
behavior is shown in the phase diagrams Fig. [T]and[3l 

This analysis is motivated by a conjecture by Susskind [3] on a correspondence between long 
fundamental strings and a small black hole, and it is thus of interest to understand the behavior of 
the size of a long string, especially when it contracts. In the case of pure Newton gravity, Khuri[T6] 
has carried out a variational analysis for Edwards Hamiltonian, and derived a scaling consistent 
with the result of Horowitz and Polchinski [5]. In this paper, a generalization with the repulsive 
force is considered, and we have estimated the critical coupling and the critical size, at which a 
string may be transformed into a corresponding black hole. The repulsive force pushes the critical 
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gravitational coupling to larger values, and the size of a black hole at the transition point will 
swell, if the repulsive force is sufficiently strong. Although the origin of the repulsive force in string 
theory is not transparent, it has been argued that it would emerge nonperturbatively to explain 
the exponential spreading of string degrees of freedom in the context of black hole complementarity 
[14t 120] , and our analysis may serve an interesting observation when the repulsive force is at 
work. There might be other sources to induce effective self-avoiding property; for example, i?-field 
(or other higher rank tensors) exchange or fermionic degenerate force. The former would not be 
important for spherically symmetric configurations as it cancels out under the average, but may 
be effective if anisotropy is introduced. The latter plays a central role in the gravitational collapse 
of stars, and if a polymer carries fermionic degrees of freedom, such as a model of superstring, the 
degenerate force can also be important to determine the properties of the polymer when the size 
is sufficiently small. 

We have been concentrating on the size of a conhguration throughout this paper. There are 
more quantities of interest such as density, or elasticity (or pressure) distribution. These are 
information necessary to write down an equation of state. Once we know these quantities, we can 
analyze the change of the size with more sophisticated methods which are used for gravitational 
collapse of stellar objects. Since a strong gravitational force is required to crush an object if a 
repulsive force is introduced, there might be a region where we need to invoke general relativity 
to describe the whole process, instead of Newton gravity. If it is the case, it will be interesting to 
investigate Tolman-Oppenheimer-Volkoff equation [21] for collapsed polymers where one can treat 
the dynamics of the system including background spacetime in a self-consistent way. 

Though our study is indeed motivated by the string/black hole correspondence, the analysis 
will also be interesting from the point of view of polymer physics. As argued in this paper, it is 
intriguing that there exist nontrivial scalings for real polymer with a long-range attractive force. 
By scaling the attractive coupling with respect to N, we can determine the scaling exponent for a 
hxed magnitude of the repulsive force {u ~ 0(1) for real polymers). For example, we take g ~ N°‘ 
and u ~ 0(1), then the size scales as i? ~ N~^ when Newton and the repulsive forces are in 
balance. It has been known that these exponents will get modified if we employ renormalization 
group analysis |12] . It will also be interesting to carry out renormalization group analysis in the 
current model to obtain more accurate exponents. 

A related observation concerns the size at the corresponding point for real polymers with 
u ~ 0(1). The size coincides with that of close-packed configurations in d dimensions, Nd, which 
is the smallest possible size scaling for A^-step SAW. In the case of real polymers, the fundamen¬ 
tal molecule structure may be broken down and the mathematical description may not be valid. 
However, it is curious that this scaling may be just a coincidence, or a more profound reason exists 
behind it. If Newton gravity is not valid before the configuration reaches this size, it is interesting 
to check whether this smallest size scaling still appears or not in analysis with general relativity. 

Another interesting issue may be about the symmetry of configurations. Although the statistical 
average leads to spherically symmetric configurations, each snapshot of configurations is known to 
be aspherical |22|. Once interaction becomes effective, the shape will be further distorted. Thus, 
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the gravitational collapse of polymers should have much richer contents than those discussed here, 
which may deserve further investigation. 

We hope to revisit these issues in future and to report elsewhere. 
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A Miscellaneous Calculations 


In this appendix, we present some details of the evaluation of each term in (12.361) and ()2.37p . We 
start with the ones involving the kinetic term, which contains a derivative. It is better to come 
back to a discrete description, da —)• Y2n=i “ Rn-i (for cr E (n — 1, n]), where 

we use a subscript to denote the position of n-th monomer. Thus, 



Rn-l)") 


/ 


I 



l)G'(n - I,n),G'(n,iV)(R„ - Kn-if 


=a^N, 


(A.l) 


and 




(R(iV) - R(0))2 


I « / \ / 

^ WR-n “ Rn-l)^ [(Rat ~ Rn) + (Rn ~ Rn-l) + Rn-l] J 

1 ^ r 

^ ^ ((Rjv - Rn)')' ((Rn - Rn-l)')' + ((Rn " Rn-l)^)' + ((Rn-l)')' ((Rn " Rn-l)')' 


n=l 

-A 1)2 


=a'^tN N + 


d ’ 


(A.2) 


where we have used R(0) = 0 and the properties of Gaussian average. 

Next, we move on to the interaction terms. As before, V{a,a') is symmetric, and we take 
cr' > cj by rewriting the integral as 2 da' da. It should be noted that we take a' > a, not 
a' > a. Since we do not consider the self-interaction of each monomer, the point a' = a should be 
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excluded. For a discrete description, we then consider n' > n. The expectation values of the delta 
function part are 


-R(u'))) = J VR{a)VR{a')VR{N)G\0,a)G'{a,a')G'{a',N)6^'^\R{a)-Yi{a’)) 

d 

d 


2TTa^G{a' — a) 


and 


(Rjv - Ro)^5^‘^^(R(ct) - R(ct')))' =a^f{N - a' + a) 


d 


(A.3) 


(A.4) 


— a) 

The Newton interaction part is again evaluated by converting to the exponential form ()2.15p 
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(A.5) 


For the one with (R(iV) — R(0))^ insertion, the expectation value is 


e-*^(^{'")-^{'"'))'(R(iV) - R(0))^ 


=a^f 


Ar / / / \ / 1 2ia^i‘^{a' — a) 

N — a + a + (a' — a) { 1 H-^- -x 


d 


-In 


2iaH^a'-a) Y^ 

^ +- d - 


(A,6) 


Upon x integral, we find 

^2 i)d—2 


/'ild-2 r^(0)) ) — v2a2/ ^ 


\R(a) - Ria' 


2-d 
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N + Y — a) 


2-d 


d 


(A.7) 


By collecting these results, we obtain (|2.38p . 


B A generic power potential and van der Waals interaction 

The methods used in this paper can be applied to any power-law long-range force. We first present 
a general formula and argue a potential problem that arises for interactions with higher inverse- 
power. As an example, we see that the size scaling due to van der Waals inverse-sextic potential 
of a real polymer in three dimensions may have difficulty. 

We consider a generic power-law potential term in d dimensions. 


U„ =■ 




R(cj) - R(cj') 


(B.l) 
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where ^ is a dimensionless coupling (positive for an attractive case), a = d — 2 corresponds to 
Coulomb-type interactions, including Newton gravity. For variational calculation, we evaluate the 
expectation value of this interaction term by use of a harmonic trial Hamiltonian, 


{Va)o = 


r(^) / qd 
r(|) \2Fi{a,a']q)J 


(B.2) 


where Fi{a,a'‘,q) is given in (I2.17p . On the other hand, in the uniform expansion model, we need 
to evaluate the following two terms with Gaussian Hamiltonian (12.3211 . 


{VaY 

(H„(R(iV)-R(0))2)' 
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N — [a — a) 
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(B.3) 

(B.4) 


In these calculations, one would have faced divergent integrals for some large values of a. This is 
due to a short distance singularities |R((t) — R((j')| —)• 0, but phenomenologically these singularities 
are avoided by using an effective repulsive force (thus a more realistic phenomenological potential 
is of Lennard-Jones type). We neglect such divergence and take only q and a dependence. 

The variational calculation for the optimized value of q leads to 


0.-2 d — 2 

1 — N^q 2 -\-Nuq 2 =0, 


(B.5) 


while the UEM consistency condition is given by 

1 - + ua-^N^ - = 0 , 


(B.6) 


where we have dropped positive numerical coefficients as before. If the repulsive force is absent 

(namely tt = 0), the variational calculation is more reliable, and the optimal value of go and the 

2 1 

scaling size are given by go — (-^0 and R ~ £(iV ^)“-2 respectively. The size increases as ^ 
gets larger if a > 2 and it is an unreasonable behavior. In the case of Newton gravity [a = d — 2), 
this corresponds to d > 4, as we have observed in Sec. 12.11 If the repulsive force is turned on and 
is balanced with Newton force, both methods lead to the same size scaling 

( 7 / \ d — OL 

• (B.7) 

Therefore, if a > d, the configuration expands if the attractive force gets stronger. In the case of 
Newton gravity, this never happens and then we obtain a reasonable scaling if the repulsive and 
Newton force balance out. However, if the inverse-power of the potential is too large, we again face 
an unreasonable behavior, and our approximations both fail. 

One of the typical example of this pathological behavior is van der Waals interaction in three 
dimensions, where a = 6 and d = 3. At this moment, it is not so clear why we fail to observe 
that van der Waals force, one of the realistic interactions of real polymers, makes a configuration 
smaller; but one possible explanation may be about the validity of the perturbative treatment of 
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the interaction, as follows. By recalling a scaling argument at the beginning of Section [2l we find 
that the marginal van der Waals coupling for a free configuration is ~ N, namely an extremely 
strong coupling, and an 0(1) coupling does not affect the configuration in three dimensions. An 
0(1) van der Waals coupling becomes effective only when the size of the configuration becomes 
R ~ . This represents close-packed configurations, and then van der Waals interaction may 

be important only for phases in dense-packing, like crystallization, for real polymers. Thus, it may 
not be justified to perturbatively treat van der Waals interaction as a long-range interaction in the 
context of determining the shape of long polymers. 

References 

[1] D. Mitchell and N. Turok, “Statistical Mechanics of Cosmic Strings,” Phys. Rev. Lett. 58 
(1987) 1577. 

D. Mitchell and N. Turok, “Statistical Properties of Cosmic Strings,” Nucl. Phys. B 294 (1987) 
1138. 

[2] J. L. Manes, “String form-factors,” JHEP 0401 (2004) 033 |hep-th/0312035|. 

J. L. Manes, “Portrait of the string as a random walk,” JHEP 0503 (2005) 070 
|hep-th/0412104|. 

[3] L. Susskind, “Some speculations about black hole entropy in string theory,” In *Teitelboim, 
C. (ed.): The black hole* 118-131 |hep-th/9309145|. 

[4] G. T. Horowitz and J. Polchinski, “A Correspondence principle for black holes and strings,” 
Phys. Rev. D 55 (1997) 6189 |hep-th/9612146 . 

[5] G. T. Horowitz and J. Polchinski, “Selfgravitating fundamental strings,” Phys. Rev. D 57 
(1998) 2557 |hep-th/9707170|. 

[6] J. L. F. Barbon and E. Rabinovici, “Touring the Hagedorn ridge,” In *Shifman, M. (ed.) et 
ah: From fields to strings, vol. 3* 1973-2008 |hep-th/0407236|. 

M. Kruczenski and A. Lawrence, “Random walks and the Hagedorn transition,” JHEP 0607 
(2006) 031 |hep-th/0508148]. 

[7] T. G. Mertens, H. Verschelde and V. 1. Zakharov, “Near-Hagedorn Thermodynamics and 
Random Walks: a General Formalism in Curved Backgrounds,” JHEP 1402 (2014) 127 
[arXiv: 1305.7443 [hep-th]]. 

“Random Walks in Rindler Spacetime and String Theory at the Tip of the Cigar,” JHEP 1403 
(2014) 086 [arXiv: 1307.3491 [hep-th]]. 

“The thermal scalar and random walks in AdS'i and BTZ,'" JHEP 1406 (2014) 156 
[arXiv: 1402.2808 [hep-th]]. 

“Near-Hagedorn Thermodynamics and Random Walks - Extensions and Examples,” JHEP 
1411 (2014) 107 [arXiv: 1408.6999 [hep-th]]. 


23 


“On the Relevance of the Thermal Scalar,” JHEP 1411 (2014) 157 [arXiv:1408.7012 [hep-th]]. 
“Perturbative String Thermodynamics near Black Hole Horizons,” arXiv:1410.8009 [hep-th], 
“The long string at the stretched horizon and the entropy of large non-extremal black holes,” 
arXiv: 1505.04025 [hep-th], 

[8] T. Damour and G. Veneziano, “Selfgravitating fundamental strings and black holes,” Nucl. 
Phys. B 568 (2000) 93 [hep-th/9907030|. 

[9] D. Amati and J. G. Russo, “Fundamental strings as black bodies,” Phys. Lett. B 454, 207 
(1999) [arXiv:hep-th/9901092|. 

[10] S. Kawamoto and T. Matsuo, “Emission spectrum of soft massless states from heavy super¬ 
string,” Phys. Rev. D 87 (2013) 12, 124001 [arXiv: 1304.7488 [hep-th]]. 

[11] Paul J. Flory, “The configuration of real polymer chains,” The Journal of Ghemical Physics 
17.3 (1949): 303-310.; 

Paul J. Flory, “Principles of polymer chemistry,” Gornell University Press, (1953) 

[12] Pierre-Gilles De Gennes, “Scaling concepts in polymer physics,” Gornell university press (1979) 

[13] Masao Doi and Samuel Frederick Edwards, “The theory of polymer dynamics,” Oxford: 
Glarendon Press, 1986. 

[14] L. Susskind, “The World as a hologram,” J. Math. Phys. 36 (1995) 6377 |hep-th/9409089|; 
L. Susskind and J. Lindesay, “An Introduction To Black Holes, Information And The String 
Theory Revolution: The Holographic Universe,” World Scientific Publishing Gompany (2004) 

[15] S. Kalyana Rama, “Size of black holes through polymer scaling,” Phys. Lett. B 424 (1998) 39 
[hep-th/9710035|. 

[16] R. R. Khuri, “Selfgravitating strings and string/black hole correspondence,” Phys. Lett. B 
470 (1999) 73 [hep-th/9910122|. 

[17] R. R. Khuri, “Black holes and strings: The Polymer link,” Mod. Phys. Lett. A 13 (1998) 1407 
gr-qc/9803095|. 

[18] S. F. Edwards, and M. Muthukumar, “The size of a polymer in random media,” The Journal 
of chemical physics 89.4, 2435-2441 (1988) 

[19] P. Pfeuty, R. M. Velasco, and P. G. De Gennes, “Conformation properties of one isolated 
polyelectrolyte chain in D dimensions,” J. Phys. (Paris) Lett. (Journal de Physique Lettres) 
38, 5-7 (1977) 

[20] K. Ropotenko, “What is the rate at which entropy of a string falling toward a black hole 
increases?,” Phys. Rev. D 79 (2009) 064003 [arXiv:0809.5236 [hep-th]]. 


24 


[21] R. C. Tolman, “Effect of imhomogeneity on cosmological models,” Proc. Nat. Acad. Sci. 20 
(1934) 169 [Gen. Rel. Grav. 29 (1997) 935]. 

R. G. Tolman, “Static solutions of Einstein’s field equations for spheres of fluid,” Rhys. Rev. 
55 (1939) 364. 

J. R. Oppenheimer and G. M. Volkoff, “On Massive neutron cores,” Phys. Rev. 55 (1939) 374. 

[22] Joseph Rudnick and Geoge Gaspari. “The aspherity of random walks,” Journal of Physics A: 
Mathematical and General 19.4 (1986): L191. 

S. J. Sciutto, “Study of the shape of random walks,” Journal of Physics A: Mathematical and 
General 27.21 (1994): 7015. 

Gharbel Haber, Sami Alom Ruiz, and Denis Wirtz, “Shape anisotropy of a single random-walk 
polymer,” Proceedings of the National Academy of Sciences of the United States of America 
97.20 (2000): 10792. 


25 



